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Abstract 

We analyze the interior geometry of static, spherically symmet- 
ric black holes of the Einstein- Yang-Mills-Higgs theory. Generi- 
cally the solutions exhibit a behaviour that may be described as 
"mass inflation", although with a remarkable difference between 
the cases with and without a Higgs field. Without Higgs field the 
YM field induces a kind of cyclic behaviour leading to repeated 
cycles of mass inflation - taking the form of violent explosions - 
interrupted by quiescent periods and subsequent approaches to an 
almost Cauchy horizon. With the Higgs field no such cycles occur. 
In addition there are non-generic families with a Schwarzschild 
resp. Reissner-Nordstr0m type singularity at r = 0. 
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1 Introduction 



Most of the knowledge on the interior geometry of static black holes derives 
from the well-known exact solutions like the Schwarzschild (S) and Reissner- 
Nordstr0m (RN) solution. Whereas for the S-solution the radial coordinate 
stays time-like inside the horizon all the way down to the central singularity, 
the RN-solution exhibits a second - inner - horizon, behind which the radial 
coordinate becomes space-like again. This inner horizon is a so-called Cauchy 
horizon, beyond which the evolution of any matter system is influenced by 
data causally disconnected from the earlier history Jl|. Furthermore, while 
the outer horizon is an infinite redshift surface, the Cauchy horizon is a 
surface of infinite blueshift. This suggests, that the Cauchy horizon should 
be unstable against small perturbations and a space-like singularity should 
form replacing the former Cauchy horizon ||. In fact, perturbative studies 
indicate an exponential growth of the mass-function close to the inner horizon 
- a phenomenon dubbed "mass inflation" |2j . However, as yet no unanimous 
conclusion about the reliability and genericity of these perturbative results 
seems to have been obtained |J. In this situation it is definitely of interest to 
investigate the internal structure of black holes with other forms of matter 
like Yang-Mills and Higgs fields. This is the subject of the present study. 
Compared to the previously mentioned perturbative analysis our work is 
"exact" , although at least partially based on numerical results. 

One may summarize these results saying that generically no Cauchy hori- 
zon forms, because its existence requires fine-tuning of the initial data at the 
outer horizon. Solutions with a Cauchy horizon are in fact obtained through 
such fine-tuning, leading to a RN-type singularity at r = 0. Generically, 
however, the solutions exhibit "mass-inflation" , although with a remarkable 
difference between the cases with and without a Higgs field. The behaviour is 
actually much simpler with the inclusion of a Higgs field. In this case the mass 
function shows generically exponential growth all the way down to r = 0. 
On the other hand in the seemingly simpler case without the Higgs field the 
situation is more involved. The YM field induces a kind of cyclic behaviour 
leading to repeated cycles of mass inflation - taking the form of violent "ex- 
plosions" - interrupted by quiescent periods and subsequent approaches to 
an almost Cauchy horizon. This behaviour is particularly spectacular due to 
the fantastic growth of the mass function during these explosions by hundreds 
of orders of magnitude, a phenomenon unprecedented in standard physical 
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problems. Actually, these explosions become exponentially more violent af- 
ter each cycle such that it is practically impossible to follow more than one 
or two of them numerically. 

A further novel phenomenon is a kind of chaotic behaviour generated 
by a hierarchical structure of families of non-abelian RN-type (NARN) and 
non-abelian S-type (NAS) solutions separating generic solutions with cyclic 
(EYM) or acyclic (EYMH) behaviour. 

Our numerical results exhibiting these claims are predominately based 
on the EYM system, but we have no doubt that we could easily establish 
similarly convincing evidence in the EYMH case. 

In this discussion the global behaviour outside the horizon was ignored. 
Clearly it is of interest to find out what is the behaviour of asymptotically flat 
solutions, which for the EYM theory constitute a discrete set of 1-parameter 
families - described by curves in the (Wh, r h ) plane ||, 15] . Obviously without 
further restriction they will show the generic behaviour inside the horizon, 
although for certain discrete points S-type resp. RN-type singularities at 
r = are possible. 

In a recent paper Donets et al. presented their results on the internal 
structure of static, spherically symmetric black holes of the Einstein- Yang- 
Mills (EYM) theory (without Higgs field). Our results for this particular 
case essentially agree with the findings of Donets et al., although we differ 
in some details. In particular the chaotic stucture we find close to the basic 
NARN solutions was not observed in ||. 



2 Ansatz and Field Equations 

The contents of this chapter are essentially (up to minor changes in notation) 
a copy from our earlier paper |7j], which we include for the convenience of the 
reader. 

For the static, spherically symmetric metric we use the parametrization 
ds 2 = A 2 Bdt 2 - ^ - r 2 (R)dQ 2 , (1) 

with dQ 2 = d6 2 + sin 2 ^^ 2 and three independent functions A, B, r of a 
radial coordinate R, which has, in contrast to r, no geometrical significance. 
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It is common to express B through the "mass function" m defined by B = 
1 — 2m/r. 

For the SU(2) Yang-Mills field we use the standard minimal spheri- 
cally symmetric (purely 'magnetic') ansatz 

W^dx^ = W(R) {T x dO + T 2 sin 9d<p) + T 3 cos 9d<p , (2) 

and for the Higgs field we assume the form 

$ a T a = H(R)n a T a , (3) 

where T a denote the generators of SU(2) in the adjoint representation. Plug- 
ging these ansatze into the EYMH action results in 

\(l + B{{r'f + ^V)')) " BrW, - V 2 ] , (4) 



s = - i (in a 

with 
and 



V^ + \t»>?, (5) 



(I _ w 2 ) 2 B 2 r 2 
V 2 = [i Jl j + ^-(H 2 - a 2 ) 2 + W 2 H 2 . (6) 
2r z 8 

Through a suitable rescaling we have achieved that the action depends only 
on the dimensionless parameters a and f3 representing the mass ratios a = 
MwVG/g = Mw/gM-pi and (3 = Mh/Mw (g denoting the gauge coupling 
and G Newtons constant). Q 

We still have to choose a suitable gauge for the radial coordinate R. The 
most natural choice is R = r, i.e. Schwarzschild (S) coordinates. However, 
this choice is singular at a stationary point of r(R) ('equator'). Such singu- 
larities are avoided using the gauge B = r~ 2 for B > resp. B = — r~ 2 for 
B < 0. We denote this radial coordinate by r in order to distinguish it from 
the S coordinate r. With this choice the spatial part of the metric takes the 
simple form 

ds 2 = r 2 (dr 2 + dn 2 ) , (7) 



The usual quantum definition of Mpi = y is obtained from our classical expression 

replacing the dimensionful classical gauge coupling g 2 by the dimensionless j- and putting 
the latter equal to one. 



4 



suggesting to call them isotropic coordinates. 

Using S coordinates the field equations obtained from (P are 



(BW')' = W{ 



W 2 



+ H 2 ) - 2rBW'V 1 



2^2 



(3 2 r 



(r 2 BH')' = (2W 2 + ~^-(H 2 - a 2 ))H - 2r 3 BH% 

(rB)' = 1 - 2r 2 BV 1 - 2V 2 , 
A' = 2rV x A . 



(8a) 

(8b) 

(8c) 
(8d) 



The equations obtained with isotropic coordinates B = — r 2 are essentially 
Eqs. (9) of [0, there are however some sign changes due to B < 0. 



r 

N' 

K 

w 

u 

H 

V 



rN , 

: (k - N)N - 2U 2 - V 2 

f3 2 r 2 



-1 - k 2 + 2U 2 



(H 2 - a 2 ) 2 + 2H 2 W 2 



■ rU 



W{W 2 - V 



rH 2 W - (« - N)U 



V 



2J2 



(3 2 r 



■{H 2 - a 2 )H -2W 2 H - kV 



(9a) 
(9b) 

(9c) 
(9d) 

(9e) 
(9f) 

(9g) 



together with the constraint 



and the definitions 

N=-, 
r 



2kN = -1 + iV 2 + 2U 2 + V 2 + 2V 2 , (10) 



(A 2 B)' W 
K =L^±-+N, U = — , and V = H'. (11) 



3 Singular points 

The field Eqs. (§) are singular at r = 0, r = oo and for points where B 
vanishes. Whereas the former singularities are of geometrical origin (the 
action of the rotation group degenerates) the zeros of B turn out to be 
coordinate singularities, in fact, of two different kinds: 
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1. Points where B = 0, but A 2 B ^ are stationary points of the function 
r(R) and hence choosing coordinate leads to singular derivatives; 
our second coordinate choice avoids this problem. 

2. Points where B = and A 2 < oo correspond to an horizon; as is well 
known a retarded time coordinate t* —t — J avoids this singularity. 

As we want to stick to the time coordinate t we have to treat horizons as 
singular points. In order to guarantee the fmiteness of A we have to require 
the regularity conditions 



1 dVo 

rB'\ h = l-2V 2 \ h , B'W'L 



2 dW 



h 



r B H \ h 



dH 



(12) 



In general it is not possible to parametrize regular solutions directly by 
their initial data at the singular point. However, as was already discussed in 
[[], 0] in the present case this is in fact possible, using a sharpened version 
(Prop. 1 of ||) of the standard text-book existence theorems 0. 
Proposition: 

Consider a system of first order differential Eqs. for m+n functions y = (u, v) 

s^ = sfi(s,y) , i=l,...,m, (13a) 

S ~T = + s 9i(s,y) , i = l,...,n, (13b) 
as 

with constants A» > and let C be an open subset of R m such that /j, are 
analytic in a neighbourhood of s — 0, y — (c, 0) for all c G C. There exists 
an m-parameter family of local solutions y c (s) analytic in c and s for cGC, 
\s\ < Sq(c) such that y c (0) = (c, 0). 

In order to meet the requirements of this Proposition we may introduce 
the coordinate s = r — r h and put M 

ui = r , u 2 = W, u 3 = H, (14a) 

ui = - - — ( 1 - 2V 2 ) , (14b) 



BW 



it' 



V 2 = «2Hl-+«3 > ( 14C ) 



^3 



I^-^tii + £«?(«!-!)). (14d) 
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There remain two parameters - Wh and Hh - to describe solutions with a 
regular horizon at r^. For an event horizon r has to increase as one deviates 
from the horizon in the direction where B is positive. This requires the 
inequality 

rB'\ h = l-2V 2 \ h > , (15) 

which for the case without Higgs field reduces to r h > \W% — 1|. For V 2 \h = 
1/2 we get B'\h = 0, which looks like the condition for a degenerate horizon. 
As discussed in |7j the latter is, however, only obtained for 1 — 2V 2 \h = 
dV 2 /dW\h = dV 2 /dH\h = 0; without Higgs field this is satisfied only for 
rh = 1, Wh = yielding the extremal Reissner-Nordstr0m solution. Initial 
data with V 2 \ h = 1/2 but dV 2 /dW\ h ^ and/or dV 2 /dH\ h ^ turn out 
to describe a non-degenerate horizon, which is simultaneously a maximum 
of r and therefore has a singularity in S coordinates. In order to avoid this 
singular behavior one may use isotropic coordinates, in which the boundary 
conditions at the horizon (assumed to be at r = 0) read 



(16a) 
(16b) 
(16c) 
(16d) 




N(r) = N 1 (r + 0(r 3 )) - ( W 2 X + ^-H\ ) r 3 + 0(r 5 ) , 



T 



W(r) = W h + r h W 1 - + 0(r 4 ) , 



T 



H(r) = H h + H 1 - + 0(T 4 ) 



where 



AT, 



2 + 



= — r 



4 avy 



^1 



1 ay 9 



2 <9# 



(17) 



The function k(t) has a simple pole at the horizon, but k(t) — 1/r is regular, 
i? = — iV 2 , PF, and if are analytic in r 2 and r 2 is analytic in the S coordinate 
r as long as V 2 \h 7^ 1/2 (i.e., #1 / 0). For the special case V 2 \h = 1/2 we 
obtain 



B 



-8(W? + F 2 /2)5(l - r /r,)5 + 0((1 - r/r,) 2 ) 



W=W fc + r fc Wi 



r/r h 



W? + H 2 /2 / 



+ 0(l-r/r h ) 



(18a) 
(18b) 
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ff=fl ' +ffl (i^)' +0(1_rw ' (18c) 

In the following we will use the terms event resp. Cauchy horizon for any 
horizon with B'\h > resp. < 0, although the original meaning of these 
terms applies only to asymptotically flat solutions. 

Next we turn to the singular behaviour at r = 0, which is of particu- 
lar relevance for the internal structure of black hole solutions. We have to 
distiguish two cases, B > and B < 0. 

1. B > 0: For black holes this case is only possible, if there is a sec- 
ond, inner horizon. The generic behaviour of the solutions in isotropic 
coordinates without the Higgs field is described by Prop. 13 of [|]. 
Translating to S coordinates and taking into account the Higgs field we 
get 



2{L — W , 



H(r)=H + H ir + 0(r 2 ) , (19b) 

fl(r)= 0M_^ + o(1) , (19c) 



r 



which is a 5-parameter, i.e. generic, family of solutions. In the case 
with vanishing Higgs field Prop. 13 of [BJ implies the analyticity of 
W(r) and r 2 B(r), i.e. the expressions above are in fact the beginning 
of a convergent Taylor series. More generally, in the case with a Higgs 
field, introducing the variables 

Ul = W , u 2 = (rBW + W(W 2 - l))/r , (20a) 
u 3 =(ff^-l)/r, (20b) 
u 4 = H , u 5 = H' , (20c) 

it is easy so check that 

ru£ = 0(r), i = l... 5, (21) 
and thus analyticity follows from the Proposition with C = {Wq ^ 1}. 
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According to the asymptotics of B(r) we may call the singular be- 
haviour to be of RN-type. The special case Wq = 1,M < leads to 
S-type behavior with a naked singularity. On the other hand Wq = 
1, H = 0, M = gives regular solutions. 

2. B < 0: This case is more involved, with two disjoint families of 
singular solutions. 

2.1 There is a 3-parameter family of solutions with a S-type singularity, 
characterized by the asymptotics 

W(r) = 1 + W 2 r 2 + 0(r 3 ) , (22a) 
H(r) = H + O(r) , (22b) 

B(r) = -— + 0(1). (22c) 



Introducing the variables 

rB 



Ul = BW' , u 2 = \, u 3 = H, (23a) 



(1 - W 2 ) w 
r 2 2r 

v 2 = rBH' - 2H , (23c) 



«i = - ^ + — , (23b) 



we obtain 



r^ = 0(r), fori = 1...3, (24a) 
rv'i = -2v! + 0(r) , rv' 2 = -v 2 + 0(r) . (24b) 

With C = {M > 0} the analyticity follows again from the Proposition. 

Obviously the condition B < (Mo > 0) prevents the existence of 
regular solutions in this case. 

2.2 There is an additional 2-parameter family of solutions with a pseudo- 
RN singularity (pseudo because B < 0). 

W(r) = W ±r + 0(r 2 ) , (25a) 
H(r) = H + O(r 2 ) , (25b) 
Bjr) = JW-l)\mil-WS) +m (25c) 



Introducing the variables 



u x = W 



we get 



Vi : 
v 2 \ 



rv'n 



1 =F [W 1 - 

(W 2 - l) 2 
r^B 



u 2 = H , 

W{W 2 - 1] 



rB 
l)/r, 



/r, u 3 = #', 



O(r) fori = 1,2 
-2^i -4w 2 + 0(r) 
-w 3 + 0(r) . 



rv' 2 



Vi — v 2 + 0(r) 



(26a) 
(26b) 

(26c) 



(27a) 
(27b) 
(27c) 



With C = {Wq 7^ 1} the Proposition implies again analyticity. The 
negative eigenvalues of the linearized equations are Ai j2 = —1/2(3 ± 
i^/l5) and A3 = — 1. This is a repulsive focal point which will turn out 
to be important for the cyclic behaviour in the EYM case. 

In the case B < we obtained no singular class that has enough parameters 
(three for EYM and five for EYMH) to describe the generic behaviour. Since, 
also the appearance of a second, inner horizon is a non-generic phenomenon, 
one may wonder, what the generic behaviour inside the horizon near r = 
looks like. This question actually is at the origin of our interest in this 
problem. Some insight was obtained by a detailed numerical study, whose 
results will be presented in the next chapter (compare also 0). 

Finally we would like to address the question of geodesic incompleteness 
inside the horizon expected in view of general singularity theorems PJ. We 
will show that for solutions without a Cauchy horizon (i.e., B < for r < r^) 
the radial time- and lightlike geodesies are incomplete. Their equations are 



fdt\- 



g 



dr\ 2 



with 



1 or . 



(28) 



From staticity we get the constancy of gudt/dr = E and thus 

/dr\ 2 E 2 



I dry _ E z 



eB . 



(29) 
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Figure 1: First inflationary cycle of the fundamental black hole of EYM 
theory with = 1 and Wh = 0.6322 . 

From Eq. (|8]d) we see that A is a monotone function of r implying 



where is the difference in proper time resp. affine parameter between the 
horizon and the origin. 

4 Numerical results 

In order to investigate the generic behaviour of non-abelian black holes inside 
the event horizon, we integrate the field Eqs. (§) resp. (|9|) from the horizon 
assuming B < 0, ignoring the constraints on the initial data at the horizon 
required for asymptotic flatness. 

As mentioned before, the horizon is a singular point of the equations. 
Consequently one has to desingularize the equations in order to be able to 
start the integration right there. How this can be done, was described in the 
previous chapter resp. in J7|, §. 

If one ignores the behaviour of the solutions outside the horizon, one may 
relax the condition (|T5| ) ensuring a decreasing resp. increasing r inside resp. 
outside the horizon. At least part of the solutions violating (pl|) possess an 
equator, i.e. a maximum of r. For those the use of S coordinates is excluded. 




(30) 
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Figure 2: First two cycles of the solution with r/ l =0.97 and Wh = 0.2. For 
the second cycle a suitably stretched coordinate x is used. 



As one performs the numerical integration one quickly runs into serious 
problems due to the occurence of a quasi-singularity, initiated by a sudden 
steep raise of W and subsequent exponential growth of B resp. m (compare 
Figs. [1], |2|, and |3| for some examples). This inflationary behaviour of the 
mass function is similar to the one observed for perturbations of the RN 
solution at the Cauchy horizon 0. While this exponential growth continues 
indefinitely for the EYMH system, it comes to a stop without the Higgs field. 
The mass function reaches a plateau and stays constant for a while until it 
starts to decrease again. When B has become small enough, i.e. the solution 
comes close to an inner horizon, the same inflationary process repeats itself. 
Generically this second "explosion" is so violent (we will give estimates on 
the increase of m in chapter |5|) that the numerical integration procedure 
breaks down. 

Besides these generic solutions there are certain families of special solu- 
tions obtained through suitable fine-tuning of the initial data at the horizon. 
There are two classes of such special solutions. The first class are black holes 
with a second, inner horizon, the second are solutions with one of the sin- 
gular behaviours at r = for B < described in chapter [| The numerical 
construction of such solutions is complicated by the fact that both boundary 
points are singular points of the equations. The strategies employed to solve 
such problems are well described in our paper on gravitating monopoles [/]]. 
Actually, in order to control the numerical uncertainties we used two differ- 
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Figure 3: Inflationary solutions with Higgs fields with a = 0.2, (3 = 0, 
r h = 1.2 and W h = 0.15, if h = 0.8 resp. W h = 0.2, H h = 0.5 . 

ent methods, which may be called "matching" and "shooting and aiming". 
For matching we integrate independently from both boundary points with 
regular initial data, tuning these data at both ends until the two branches 
of the solution match. For shooting and aiming we integrate only from one 
end and try to suppress the singular part of the solution at the other end by 
suitably tuning the initial data at the starting point. 

As already said, the first class of special solutions consists of black holes 
with a second, inner horizon; let us call them non-abelian RN-type (NARN) 
solutions. As was explained in chapter ^ a regular horizon requires the si- 
multaneous vanishing of B, BW and BH'. We have determined two such 
1-parameter families for the EYM system, shown in Fig. f|, whose signifi- 
cance will be explained in chapter |5|. As may be inferred from Fig. [|, the 
(dotted) curve 2 corresponding to one such family intersects all (solid) curves 
describing asymptotically flat solutions except the one for n = 1. In contrast 
to what is claimed by Donets et al. our curve continues straight through 
the parabola r h = 1 — W% and runs all the way to r h = Wh = 0. As already 
mentioned the branch to the left of the parabola cannot be obtained using 
S coordinates. The corresponding curve of Donets et al. makes a suspiciously 
sharp turn very close to the parabola and runs to the point = 1, Wh = 0. 
We are convinced that the latter piece of their curve is an artefact of unre- 
liable numerics caused by the use of S coordinates becoming singular at the 
parabola (compare the discussion in chapter ^). 
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0.5 1 1.5 2 r 2.5 

1 h 



Figure 4: Initial data for special solutions. The solid curves represent 
asymptotically flat solutions with n zeros of W. The other curves represent 
various NARN and NAS families. 

As already stressed B' has to be positive at an event horizon. This condi- 
tion is violated for values of < r p , where r p denotes the value of where 
the NARN curve 2 intersects the parabola. However it turns out to be ful- 
filled for the second horizon if < r m with r m rs 0.112 on the NARN curve. 
Thus this curve may be divided into three pieces according to < < r m , 
r m < r h < r p and r p < r^. On the first interval the two horizons have ex- 
changed their roles, whereas on the second interval both of them are Cauchy 
horizons with a maximum of r (equator) in between. The only difference 
between the first and the third interval is that Wh > 1 resp. < 1 on the 
event horizon. 

Our second NARN family (curve 5 of Fig. ^) stays completely to the left 
of the parabola and ends at w 0.9 close to the curve 3, whose significance 
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Figure 5: A NARN solution and two accompanying NAS solutions with 
W h = 0.3 and r h = 0.198728 (on curve 5 of Fig. §, r h = 0.290837 (on 
curve 6 with W — ► +1), and = 0.170525 (on curve 7 with W — > —1). 

will be explained below. 

The second class are solutions without a second horizon (i.e. B stays 
negative) approaching the center r = with one of the two singular be- 
haviours described in chapter |3], i.e. those with a S-type singularity resp. with 
a pseudo-RN-type singularity; let us denote them NAS resp. NAPRN solu- 
tions. We have determined several NAS families represented by the dashed- 
dotted curves of Fig. f|. The curve 1 staying to the right of the parabola 
coincides with the corresponding one of Donets et al., whereas the others, 
staying essentially to the left of the parabola are new. As will be explained 
in chapter |], the two NAS curves 6 and 7 accompanying the (dotted) NARN 
curve 5 are expected to merge with the NAS curve 3 close to = 0.9. Some 
of the NAS curves (e.g., 3 and 4) are expected to extend indefinitely to the 
right, but numerical difficulties (too violent "explosions") prevented us from 
continuing them further to larger values of r^. They will intersect the (solid) 
curves for asymptotically flat solutions with n = 2,3,... zeros of W and 
therefore yield additional asymptotically flat NAS black holes beside those 
found by Donets et al. ||, contradicting their uniqueness claim. 

Finally there are the NAPRN solutions, which constitute a discrete set 
according to the number of available free parameters at r = 0. We found 
several such solutions with ly'(O) = — 1 (compare Tab. p. Only one of them 
has no maximum of r and was also found by Donets et al. 
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Table 1: Initial data for several NAPRN solutions with W(0) 





Wo 




w h 


# of zeros 


1 


0.9663634 


0.7867834 


1.54085646 





2 


0.93306559 


1.889087974 


0.29197873 





3 


0.9215847 


0.8032373 


-2.06670276 


1 


4 


0.9150108 


0.009594267 


-0.8037193 


1 



5 Qualitative Discussion 



We shall now give a qualitative picture of the solutions and try to explain 
our numerical results. Since the generic behaviour of the solutions is rather 
different in the cases with and without Higgs field, we shall treat the two 
cases separately. Let us first concentrate on the case without Higgs field. 
For simplicity we introduce the notation U = BW and B = rB and use 
a = — ln(r) as a radial coordinate. Observe that B w —2m for small r. 
With these variables the field Eqs. (H) become (a dot denoting d/da) 



W 

u 

B 



2U 
-"B 



-W- 



w 2 - 1 



W 2 ) 2 



2U 3 

-l Or 

f B* ' 

- 1 ) + 2r 2 



B 



(31a) 
(31b) 
(31c) 



Close to the horizon the first term in the equation for B dominates (since 
U vanishes at r = r^) and thus B becomes negative. Provided W 2 does not 
tend to 1, this term will, however, change sign as r decreases and B will turn 
back to zero. Assuming further that U does not tend to zero simultaneously, 
the second term in the equation for U will grow very rapidly as B approaches 
zero, leading to a rapid increase of U . This in turn induces a rapid growth of 
B (compare Fig. [l]). Once the second terms in Eqs. ([3l|b,c) dominate one gets 
(U / B) and thus U/B = W'/r tends to a constant c. As long as (rc) 2 
is sizable U and B increase exponentially, giving rise to the phenomenon 
of mass inflation. Eventually this growth comes to a stop when (rc) 2 has 
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become small enough. Then U and B stay constant until the first terms 
in Eqs. (j3l[b,c) become sizable again. As before B tends to zero inducing 
another "explosion" resp. cycle of mass inflation (compare Fig. [2]). 

In the discussion above we made two provisions - that W 2 stays away 
from 1 and that U does not tend zo zero simultaneously with B. If the 
first condition is violated, i.e. W 2 1 we get a NAS solution. If on the 
other hand U and B develop a common zero we get a NARN solution, i.e. a 
solution with a second horizon. Both these phenomena can occur after any 
finite number of cycles, giving rise to several NAS resp. NARN curves as in 
Fig. f|. Generically W changes very little during an inflationary cycle, with 
the exception of solutions that come very close to a second horizon, i.e. close 
to a NARN solution. In this case W may change by any amount, depending 
on how small U becomes at the start of the explosion. By suitably fine- 
tuning the initial data at the horizon one can then obtain new NAS solutions 
with W — > ±1 or a new NARN solution. In this way each NARN solution 
is the 'parent' of two NAS and one NARN solution. Fig. f| shows two such 
generations: the NARN solutions labelled 2 have the NAS children 3 and 4 
and the NARN child 5; the curves labelled 6 and 7 are the NAS children of 5 
(see Fig. pi). Whenever the value of W at the second horizon of a NARN 
solution approaches ±1 this NARN curve and its NAS children merge with 
the corresponding sibling NAS curve having one cycle less. This hierarchy 
of special solutions gives rise to a kind of chaotic structure in this region of 
'phase space" . 

After this qualitative discussion we would like to present a simplified 
quantitative discussion of the solutions, following essentially one complete 
cycle. This will also provide us with a discrete map of the variables W, U 
and B from one plateau to the next. Let us denote their initial data at 
some plateau by Wq, Uq and B . Since the plateau is characterized by the 



(effective) vanishing of the second terms in Eqs. (31b, c) involving the factor 



r and the constancy of W, we can easily integrate them to 

u = Uq _^o}^o tl : and B = Bq+ ] , (32) 

ignoring the term — r in the B equation. The beginning of the subsequent 
explosion is characterized by B « 0, i.e. r « tq = — (Wq—1) 2 /Bq. Depending 
on the value of Wq the value of U at this point is essentially given by Uq or 
—Wo(Wq — l)/r . Whatever it is, let us denote this value again by Uq. For 
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the description of the solution through the explosion it is sufficient to solve 
the simplified system 

u u 3 u 2 

W = -r 2 ^, ^ = 2r 2 ^, B = 2r^, (33) 

implying U/B = c with some constant c. Plugging c back into Eqs. (|33|) we 
can integrate them to get 

U = U exp(c 2 (r 2 -r 2 )^) , (34a) 
B = ^exp(c 2 (r 2 -r 2 )) , (34b) 



c 



W = W + ^(r 2 -r 2 ) . (34c) 

We still have to determine c joining this solution to the one before the ex- 
plosion. Equating the derivatives of B at r = r using the expressions in 



Eqs. (H) and fl3|) we obtain c = { 2 ° Uor s 

In order to obtain the values of Wx,Ux, B\ after the explosion we may 
savely put r = in Eqs. (0) and obtain 

Ux = t/ e (cro)2 , B x = ^e (cro)2 , (35a) 

c 



c 



(w 2 -if n (w 2 - r 



2 



W = W - -r 2 , with r = — - , c = v / . (35b) 



5 2C/ r, 



o'o 



It is instructive to illustrate these relations on an example. We take the 
fundamental black hole solution with = 1 and Wh = 0.6322 shown in 
Fig. |l| For the first explosion one finds the parameters ro ~ 2.7 ■ 10~ 4 and 
c « 1.1 • 10 5 yielding cr « 30 and thus Bx ~ e 900 and W\ — W ~ 4 • 10~ 3 . 
The subsequent explosion will then take place at the fantastically small value 
r ~ e" 900 w 10" 330 . 

Since the change of W in one inflationary cycle has an extra factor ro 
the function W stays practically constant. If we furthermore concentrate on 
cases, where the first term in Eq. (j3T[b) can be neglected we may use the 
simplified system 

r7 3 • d - w 2 ) 2 u 2 

W=0, U = 2r 2 ", B= [l W) +2r 2l L, (36) 
B z r B 
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also discussed by Donets et al. 0. Introducing the variables x = rU / B = W 
and y = —(1 — W 2 ) 2 jrB one obtains the autonomous system 

W = 0, x = (y-l)x, y = y(y + l-2x 2 ). (37) 

Since the first of these equations may be ignored, we can concentrate on the 
x, y part. As usual for 2-dimensional dynamical systems the global behaviour 
of the solutions can be analyzed determining its fixed points. Since the "large 
time" behaviour a —>■ oo corresponds to the limit r — > these fixed points 
are related to the singular solutions at r = discussed in chapter |. There 
are essentially three different fixed points. 

1. For y < there is the fixed point x = 0, y = — 1 giving the RN 
type singularity. Its eigenvalues are —1 and —2, hence it acts as an 
attracting center for a —>■ oo. 

2. Then there is the point x — y = 0, a saddle with eigenvalues ±1. 

3. In addition there are the points x = ±1, y = 1 with the eigenvalues 
1/2(1 ± «a/T5), related to the pseudo-RN type singularity. This fixed 
point acts as a repulsive focal point, from which the trajectories spi- 
ral outwards. Since solutions of the approximate system given by the 
Eqs. (|3"7D cannot cross the coordinate axes, solutions in the quadrants 
y > 0, x > resp. x < stay there performing larger and larger turns 
around the focal point coming closer and closer to the saddle point 
x = y = without ever meeting it. As observed by Donets et al. this 
nicely explains the cyclic inflationary behaviour of the solutions in the 
generic case. 

Finally we come to the black holes with Higgs field. Apart from the 
generic solutions there are the special ones approaching r = with a singular 
behaviour described in chapter |3[ On the other hand, the generic behaviour 
is much simpler than in the previously discussed situation without Higgs 
field. An easy way to understand this difference is to derive the analogue of 
the simplified system Eqs. fl3"T|). Introducing the additional variable z = —H 
and ignoring again irrelevant terms one finds 

W = 0, H = -z, (38a) 
x = (y — l)x, z = yz (38b) 
y = y(y + l-2x 2 -z 2 ) . (38c) 
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Leaving aside the decoupled equations for W and H one may study the fixed 
points of the (x, y, z) system. For z = one clearly finds the previous fixed 
points of the (x, y) system. However, for z 7^ the focal point disappears 
and the only fixed point for y > is x = y = 0, z = z Q with some constant 
Zq. For Zq < 1 this point is a saddle with one unstable mode, whereas for 
Zq > 1 it is a stable attractor. The latter describes the simple inflationary 
behaviour described in chapter [| and shown in the left part of Fig. |3|. So- 
lutions approaching a fixed point with z\ < 1 eventually run away from it 
again and ultimately tend to one with z\ > 1 as shown in the right part of 
Fig. |3|. In both cases the Higgs field grows logarithmically as r — > 0. 

Finally let us remark that all these solutions are singular at r = 0, since 
a regular origin requires B(0) = 1. 



6 Summary 

We study the behaviour of the static, spherically symmetric non-abelian black 
holes inside their event-horizon. We have chosen a Higgs field in the adjoint 
representation of SU(2), but we expect similar results for Higgs fields in other 
representations, because for r — > the differences do not seem to be relevant. 

In the generic case the solutions inside the horizon show a phenomenon 
known as "mass inflation" from perturbations of the Reissner-Nordstr0m 
black holes at their Cauchy horizon. Whereas for the black holes with a Higgs 
field the mass function grows monotonously once the inflation has started 
those without Higgs field run through an infinite sequence of inflationary 
cycles, which take the form of ever more violent "explosions". Besides the 
generic solutions there are exceptional ones with a Schwarzschild or Reissner- 
Nordstr0m type singularity at the center of symmetry. These exceptional 
solutions accumulate in a way that may be interpreted as a kind of "chaotic" 
behaviour. 

Since the typical length scale of these black holes is the Planck length (at 
least for values of the gauge coupling g 2 /Tic of order one) one may question 
the physical relevance of these results of the classical theory. However, there 



are claims |L0| that the phenomenon of mass inflation connected with the 
Reissner-Nordstr0m black hole persists in the quantized theory. Definitely 
this subject requires further study. 

Apart from that, our results provide a non-perturbative confirmation of 
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the mass inflation phenomenon observed for perturbations of the RN black 
hole. Furthermore we believe that our results are of some interest in view of 
the classical singularity theorems of Penrose and Hawking |T| . 
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